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ABSTRACT 


The  contributions  of  multiple  scattering  to  backscattered  de¬ 
polarization  ratios  in  stratus  clouds  were  investigated  using  a 
linearly-polarized  ruby  lidar  system  from  a  site  on  the  southern 
California  coast.  The  effects  of  the  receiver  field  of  view,  lidar 
elevation  angle,  and  transmitter/receiver  beam  misalignment  on 
depolarization  ratios  were  investigated. 

Results  for  a  1  mrad  transmitter  beam  width  indicate  that  a 
receiver  field  of  view  of  3  mrad  returns  an  almost  threefold  increase 
in  backscattered  depolarization  ratios  as  opposed  to  a  1  mrad  field 
of  view.  The  effects  caused  by  changing  the  elevation  angle  appear 
to  be  related  to  the  cloud  microphysical  processes  occurring  above 
the  cloud  base.  Deliberate  beam  misalignment  produces  significant 
alterations  in  the  lidar  returns,  demonstrating  the  importance  of 
maintaining  proper  alignment. 

In  the  future,  the  lidar  should  be  used  in  conjunction  with  in 
situ  microphysical  observations  and  a  program  of  theoretical  studies 
to  assess  the  degree  to  which  the  lidar  data  can  be  related  to 
specific  stratus  cloud  properties. 


TABLE  OF  CONTENTS 


Page 

ABSTRACT .  iv 

LIST  OF  TABLES . viii 

LIST  OF  FIGURES .  ix 

ACKNOWLEDGMENTS  .  xii 

Chapter 

1.  INTRODUCTION  .  1 

1.1  Background .  1 

1.2  The  Applicability  of  Lidar  .  2 

1.3  The  Stratus  Problem .  3 

2.  LIDAR  PRINCIPLES  AND  SCATTERING  THEORY  .  5 

2.1  The  Lidar  Equation .  5 

2.1.1  System  Constants  .  6 

2.1.2  The  Cross  Over  Function .  6 

2.1.3  Single  Scattering  Coefficients  .  7 

2.1.4  Multiple  Scattering  Correction 

Factor  .  9 

2.2  Depolarization  in  the  Backscatter  .  9 

2.3  Multiple  Scattering  in  Water  Clouds  .  13 

2.3.1  Early  Laboratory  and  Theoretical 

Studies .  13 

2.3.2  Early  Lidar  Measurements  .  14 

2.3.3  More  Recent  Studies .  14 

3.  STRATUS  CLOUD  PROPERTIES  .  19 

3.1  Cloud  Condensation  Nuclei  .  19 

3.2  Cloud  Water  Content  .  24 

4.  THE  FIELD  EXPERIMENT .  25 

4.1  Location .  25 

4.2  Equipment .  25 

4.2.1  Lidar  System .  25 


4.2.2  Lidar  Data  Collection  and  Analysis 
4.3  Synoptic  Situation  . 

4.3.1  Daily  Weather  Discussion  .... 

4.3.2  Satellite  Imagery  . 

4.3.3  Atmospheric  Soundings  . 

5.  EXPERIMENTAL  RESULTS  . 

5.1  Temporal  Variations  During  Stratus 

Breakup  . 

5.1.1  Visual  Cloud  Observations  .... 

5.1.2  Lidar  Findings  . 

5.2  Receiver  Beam  Width  Effects  ...... 

5.2.1  Observations  at  a  90°  Elevation 

Angle . .  •  • 

5.2.2  Observations  at  a  60°  Elevation 

Angle  . 

5.2.3  Observations  at  50°  and  30° 

Elevation  Angles  . 

5.3  Elevation  Angle  Effects  . 

5.4  Beam  Misalignment  Effects  . 

6.  CONCLUSIONS  . 

REFERENCES  . 


27 

27 

29 

29 

32 

34 


34 

36 

36 

38 

41 

45 

45 

50 

54 

58 

60 


VITA 


62 


LIST  OF  TABLES 


ap)e  Page 

2.1  Typical  microphysical  and  scattering  parameters 

for  various  water  clouds  .  8 

1  Ruby  lidar  specifications  .  28 


LIST  OF  FIGURES 


Figure  Page 


3.1  Cloud  particle  distributions  for  a  gradual  ascent 
through  shallow  nonprecipitating  maritime  stratus. 

Z  represents  height  above  cloud  base  (after  Ryan 

et  al.,  1972) .  21 

3.2  Vertical  profiles  from  stratus  clouds  advected 
onshore  during  four  periods.  Heights  relative 
to  the  fog  top  of  1.0  and  0  represent  cloud  top 
and  base  heights,  respectively  (after  Goodman, 

1977) .  22 

4.1  Map  of  southern  California  showing  the  lidar  field 

site  and  local  rawinsonde  release  locations  .  26 

4.2  The  NWS  1200  GMT  surface  chart  for  10  June  1982  .  30 

4.3  Visible  satellite  imagery  for  10  June  1982  at 

1945  GMT .  31 

4.4  Atmospheric  sounding  data  for  the  morning  of  10 
June  1982,  from  Montgomery  Field.  The  height 
above  ground  level  (AGL)  for  the  remote  sensor 

site  (185  m  MSL)  is  shown .  33 

5.1  Height  versus  time  display  of  lidar  6  values 
during  the  dissipation  of  a  stratus  layer  at 
the  indicated  times  on  10  June  1982  .  35 


5.2  Linear  depolarization  ratio  (inland  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  height  above  ground 
level  for  two  consecutive  lidar  pulses  using  1 
and  3  mrad  (MR)  receiver  beam  widths.  The  lidar 
elevation  angle  was  90°.  The  data  were  collected 

at  1036  LT  on  14  June  1982  .  40 

5.3  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 


1 


Figure  Page 


widths.  The  data  were  collected  at  a  90°  elevation 

angle  at  1234  LT  on  14  June  1982  from  a  cloud  with 

a  480  m  cloud  base .  41 


5.4  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  90°  eleva¬ 
tion  angle  at  1259  LT  on  14  June  1982  from  a  cloud 

with  a  480  m  cloud  base .  42 

5.5  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  2  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  90°  eleva¬ 
tion  angle  at  1355  LT  on  14  June  1982  from  a  cloud 

with  a  510  m  cloud  base .  43 

5.6  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  60°  eleva¬ 
tion  angle  at  1124  LT  on  14  June  1982  from  a  cloud 

with  a  450  m  cloud  base .  44 

5.7  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  60°  eleva¬ 
tion  angle  at  1217  LT  on  14  June  1982  from  a  cloud 

with  a  410  m  cloud  base .  45 

5.8  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  60°  eleva¬ 
tion  angle  at  1448  LT  on  14  June  1982  from  a  cloud 

with  a  630  m  cloud  base .  46 


x 


5.9  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  50°  eleva¬ 
tion  angle  at  1511  LT  on  14  June  1982  from  a  cloud 

with  a  580  m  cloud  base .  47 

5.10  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam 
widths.  The  data  were  collected  at  a  30°  eleva¬ 
tion  angle  at  1418  LT  on  14  June  1982.  The 

penetration  depth  scale  begins  at  540  m  AGL  .  48 

5.11  Linear  depolarization  ratio  (in  %}  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  3  mrad  observations 
taken  2  min  apart  at  the  indicated  elevation 

angles.  Average  cloud  base  was  at  545  m  AGL .  52 

5.12  Linear  depolarization  ratio  (in  %)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  vertical  penetration 
depth  above  cloud  base  for  3  mrad  observations 
taken  2  min  apart  at  the  indicated  elevation 
angles.  Average  cloud  base  was  at  545  m  AGL. 

The  data  were  collected  one  minute  later  at 

each  elevation  angle  than  in  Figure  5.11 .  53 

5.13  Linear  depolarization  ratio  (in  X)  and  relative 
returned  energy  (in  arbitrary  units)  profiles 
plotted  as  a  function  of  slant  path  range  into 
the  cloud  for  3  mrad  observations  taken  2  min 
apart  at  the  indicated  elevation  angles.  Average 
cloud  base  was  at  545  m  AGL.  The  data  were  col¬ 
lected  one  minute  later  at  each  elevation  angle 

than  in  Figure  5.11 .  55 

5.14  Effects  of  transmitter/receiver  beam  misalignment 

on  returns  from  a  stratus  cloud  layer  with  a  cloud 
base  at  320  m  above  ground  level.  Numbers  in  the 
keys  refer  to  receiver  pointing  angle  misalignment, 
in  mrad.  The  elevation  angle  was  45°  . 


56 


ACKNOWLEDGMENTS 


I  wish  to  express  my  appreciation  to  Dr.  Kenneth  Sassen,  whose 
assistance  and  guidance  made  this  thesis  possible.  I  also  wish  to 
thank  my  other  committee  members.  Dr.  Xuo-Nan  Liou  and  Dr.  Norihiko 
Fukuta  for  their  aid  and  encouragement. 

I  would  also  like  to  thank  Michael  Griffin,  Jeff  Koenig,  and 
Glenn  Ahrens  who  took  time  from  their  busy  schedules  to  assist  me 
when  the  going  was  tough. 

The  research  reported  in  this  thesis  was  supported  by  the 
National  Science  Foundation  under  Grant  ATM-82-10709.  The  lidar  data 
were  originally  collected  in  connection  with  Contract  No.  F19628-E1- 
K-0042  from  the  Department  of  the  Air  Force. 


CHAPTER  1 


INTRODUCTION 
1.1  Background 

In  1803  an  Englishman,  Luke  Howard,  first  classified  clouds 
according  to  type  by  using  their  shape,  shading,  and  general  appear¬ 
ance.  To  this  day,  a  need  still  exists  to  determine  the  exact 
composition  and  structure  of  the  various  clouds  he  classified  back 
in  1803.  Information  pertaining  to  the  amount  and  phase  of  hydro¬ 
meteors,  and  to  their  type  and  location,  can  improve  our  under¬ 
standing  of  cloud  structure  and  precipitation  development.  These 
measurements  are  invaluable  in  forming  more  precise  climatological 
and  forecasting  models  of  the  atmosphere. 

Recently,  researchers  at  the  University  of  Utah  (Sassen  et  al . , 
1985)  discovered  supercooled  liquid  water  in  cirrus  clouds,  where  in 
the  past  it  was  thought  that  only  ice  occurred.  Discoveries  such  as 
this  dictate  a  change  in  the  way  the  microstructure  and  microphysical 
processes  occurring  inside  a  cloud  are  viewed. 

Instrumentation  for  looking  at  cloud  composition  has  been 
developed  and  falls  into  three  broad  categories:  1)  in  situ  sampling 
from  airborne  platforms  (e.g. ,  aircraft,  balloons,  etc.);  2)  in  situ 
sampling  from  surface  stations  (e.g.,  rain  gauges,  mountain-top  obser 
vations,  etc.);  and  3)  remote  sensing  (e.g.,  radars,  radiometers, 
laser  radars,  etc.).  Each  system  has  unique  advantages  and 


disadvantages.  Each  specific  piece  of  equipment  is  different  and  no 
single  remote  sensor  can  supply  all  of  the  information  required  for 
atmospheric  studies. 


In  the  field  of  remote  sensing,  microwave  radars  are  the  dominant 
instruments  used.  For  many  years,  large  wavelength  radars  have  proven 
very  useful  for  looking  at  storm  structure,  and  new  techniques  have 
been  developed  for  Doppler  radars  to  infer  cloud  composition  from 
particle  fall  velocities.  However,  in  many  respects  radar  does  not 
show  as  great  a  sensitivity  to  cloud  composition  as  the  laser  radar 
(lidar) . 


1.2  The  Applicability  of  Lidar 
Over  the  past  10-15  years  there  has  been  an  increasing  use  of 
lidar  in  cloud  physics  research.  The  principle  of  lidar  remote 
sensing  is  analogous  to  radar  remote  sensing,  with  energy  radiated  in 
the  visible  spectrum  as  opposed  to  energy  in  the  microwave  region. 
Much  of  the  traditional  phemonena  associated  with  microwave  radars 
(e.g.,  the  bright  band)  exist  with  lidars  (Sassen,  1976a). 

When  a  linearly  polarized  pulse  of  high  intensity  light  is  di¬ 
rected  toward  a  cloud,  both  the  strength  and  the  linear  depolarization 
of  the  returned  signal  provide  information  on  the  phase,  type,  and 
number  or  size  of  cloud  particles.  Even  though  this  ability  exists, 
the  amount  of  lidar  research  which  has  been  published  is  not  great, 
perhaps  due  to  a  lack  of  confidence  in  the  interpretation  of  lidar 
field  data.  A  part  of  this  lack  of  confidence  arises  from  the  fact 
that  not  enough  is  known  about  the  effects  of  multiple  scattering  on 
lidar  returns.  Multiple  scattering  processes  have  been  studied  both 
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theoretically  and  experimentally  in  the  laboratory  and  field,  and 
there  is  a  large  variance  in  reported  depolarization  ratios.  It  is 
even  uncertain  to  what  extent  the  laboratory  environment  can  actually 
mimic  the  real  atmosphere.  There  is  a  need  to  better  understand  the 
physical  properties  of  the  atmospheric  target,  and  how  the  scattering 
properties  are  affected  by  the  design  of  the  lidar  system.  The  only 
way  to  do  this  is  to  use  the  lidar  in  the  atmosphere  and  compare  the 
results  to  laboratory  and  theoretical  findings.  Recently,  a  consider¬ 
able  amount  of  research  has  been  accomplished  in  the  laboratory,  and 
it  is  logical  that  the  lidar  will  again  be  used  in  the  real  atmo¬ 
sphere. 


1.3  The  Stratus  Problem 

In  the  past,  most  efforts  in  the  field  of  cloud  physics  concen¬ 
trated  mainly  on  devising  models  of  the  structure  and  development  of 
convective  clouds.  Although  convective  clouds  are  of  great  signifi¬ 
cance,  they  are  observed  much  less  frequently  than  other  cloud  types, 
such  as  stratus  and  stratocumulus.  During  the  past  20  to  30  years 
there  has  been  considerable  success  in  studying  the  dynamics  of 
stratiform  cloud  masses.  Fundamental  laws  describing  their  formation, 
development,  and  dissipation  have  been  formulated.  Satellite  obser¬ 
vations  have  also  become  increasingly  available  and  sophisticated  in 
the  past  10-15  years,  providing  valuable  descriptions  of  stratus 
clouds  on  the  synoptic  and  mesoscales.  A  better  understanding  of 
cloud  composition  and  structure  in  the  microscale  is  needed  to  go 
along  with  the  synoptic  approach. 

Stratiform  clouds  are  perhaps  most  important  when  the  radiation 
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and  energy  balance  of  the  atmosphere/ocean  system  are  considered.  As 

Randal  et  al.  (1984)  pointed  out: 

.  .  .  marine  stratus  and  stratocumulus  cloud  systems  are 
likely  to  strongly  influence  the  climate  and  climate  change. 

Their  high  albedos  (30-40%)  compared  with  the  ocean  back¬ 
ground  (10%)  give  rise  to  large  deficits  in  the  absorbed 
solar  radiative  flux  at  the  top  of  the  atmosphere,  while 
their  low  altitude  prevents  significant  compensation  in 
thermal  emission.  A  mere  4%  increase  in  the  area  of  the 
globe  covered  by  low  level  stratus  clouds  would  be  suffi¬ 
cient  to  offset  the  2-3°K  predicted  rise  in  global  temper¬ 
ature  due  to  a  doubling  of  CC^- 

The  effect  of  humanity  on  the  environment  has  been  in  the  fore¬ 
front  of  news  for  many  years.  It  is  debatable  whether  the  earth  is 
cooling  or  warming  and  what  the  effects  are  of  acid  rain,  CO2  doub¬ 
ling,  and  deforestation.  These  are  difficult  questions  and  a  better 
understanding  of  cloud  processes  is  needed  to  answer  them. 

In  conjunction  with  a  better  understanding  of  cloud  processes,  a 
better  knowledge  of  water  cloud  microphysical  properties  is  needed. 
Lidar  is  one  of  the  tools  capable  of  easily  and  economically  probing 
water  clouds.  This  research  was  undertaken  in  order  to  better  under¬ 
stand  the  multiple  scattering  process  in  water  clouds  and  particularly 
in  stratus  clouds.  The  signature  of  multiple  scattering  is  present 
in  the  backscattered  depolarization  of  the  laser  pulse,  and  it  is 
clearly  necessary  to  better  understand  this  scattering  process  in 
order  to  improve  our  knowledge  of  the  structure  of  stratus  clouds  when 
probed  with  lidar. 


As  previously  mentioned,  both  the  radar  and  lidar  have  similar 
characteristics,  including  their  governing  equations.  Battan  (1973) 
gave  a  complete  derivation  of  the  monostatic  radar  equation  and 
Sassen  (1976b)  gave  the  equivalent  derivation  for  lidar.  The  final 
form  of  the  equation  used  here  is  from  Sassen  (1977a),  where  P(R)  is 
the  power  returned  to  the  receiver  in  the  polarization  planes  both 
parallel  and  orthogonal  to  the  transmitter  as  a  function  of  range  R. 
This  equation  assumes  that  the  laser  beam  is  filled  with  targets, 
that  no  energy  scattered  out  of  the  beam  is  collected  by  the  receiver 
and  that  the  transmitter  and  receiver  beam  widths  are  equal  and  in 
parallel  alignment. 

P(R)„  ,  =~Wlf(R)  6'  (R)n  .  exp  C-2Rno(R)]  , 

11,1  8ttR^  11,1 

where  Pt  is  the  transmitted  (linearly  polarized)  power  (W),  h  the 
laser  pulse  length  (m),  Ar  the  effective  collecting  area  of  the 
receiver  (m  ),  f(R)  the  lidar  cross  over  function,  B'(R)  the  volume 
backscatter  coefficient  (km”*),  n  the  multiple  scattering  correction 
factor,  and  o(R)  the  extinction  coefficient  (km’*). 


2.1.1  System  Constants 

With  reference  to  the  seven  terms  on  the  right-hand  side  of 
Eq.  (2.1),  the  first  three  terms  in  the  numerator  are  commonly  re¬ 
ferred  to  as  system  constants.  Although  the  first  two  terms  may 
change  slightly  from  laser  pulse  to  pulse  and  so  must  be  monitored, 
the  third  term,  the  effective  collecting  area  of  the  receiver,  is 
fixed  by  the  design  of  the  lidar  receiver  (i.e.,  the  telescope). 

2.1.2  The  Cross  Over  Function 

The  fourth  term  in  Eq.  (2.1)  is  the  cross  over  function,  the 
first  of  the  range-dependent  terms.  Due  to  the  design  constraints 
of  a  high-power  monostatic  lidar  system  imposed  by  the  size  of  the 
laser,  it  is  usually  necessary  to  physically  separate  the  transmitter 
and  receiver  on  an  optical  table  and  align  their  axes  in  parallel. 
Consequently,  the  transmitter  pulse  and  the  receiver  field  of  view 
do  not  overlap  until  they  are  some  distance  from  the  lidar.  This 
distance  is  proportional  to  the  physical  separation  between  the  trans¬ 
mitter  and  receiver  on  the  optical  table,  and  to  the  laser  beam  trans¬ 
mitter  and  receiver  beam  width  angles.  Measurements  obtained  from 
targets  within  the  range  where  there  is  little  or  no  overlap  (gener¬ 
ally  a  few  hundred  meters  or  less)  must  be  examined  with  caution, 
since  the  energy  returned  from  such  short  ranges  can  only  result  from 
multiple  scattered  radiation,  assuming  negligible  side  lobe  radiation 
is  present. 


However,  in  these  experiments,  a  variable  diaphragm  was  used  so 
that  the  receiver  field  of  view  could  be  changed  between  1  and  3  mrad 


This  capability  is  a  useful  adjunct  for  multiple  scattering  studies 
since,  in  dealing  with  a  dense  target,  the  amount  of  backscattered 
power  received  with  a  lidar  becomes  a  function  of  the  beam  widths 
used.  In  effect,  this  capability  changes  f(R)  and  can  lead  to  a 
violation  of  the  lidar  equation,  Eq.  (2.1),  when  equal  transmitter 
and  receiver  beam  widths  are  not  used. 


2.1.3  Single  Scattering  Coefficients 

Since  the  first  four  terms  in  Eq.  (2.1)  can  be  considered  to  be 
constants,  the  returned  laser  power  from  atmospheric  targets  depends 
on  the  range  and  the  scattering  coefficients,  B'(R)  and  o(R).  Both 
scattering  coefficients  are  directly  proportional  to  the  cross- 
sectional  areas  and  the  particle  concentrations  per  unit  volume. 
Experimental  determinations  of  the  scattering  coefficients  have  been 
obtained,  and  typical  values  are  presented  in  Table  2.1. 

One  of  the  major  problems  associated  with  lidar  remote  sensing 
data  interpretation  has  to  do  with  deriving  the  values  of  both  B'(R) 
and  a(R).  For  radar,  under  typical  conditions,  microwave  attenuation 
is  negligible  and  the  exponential  term  in  Eq.  (2.1)  can  be  ignored. 
Radar  backscattering  is  normally  proportional  to  the  sixth  power  of 
the  diameter  of  the  particles  present,  favoring  the  detection  of  large 
precipitation  particles  within  a  cloud.  On  the  other  hand,  laser 
backscattering  and  attenuation  are  both  proportional  to  the  square  of 
the  diameter  of  the  particles  present,  leading  to  a  greater  sensi¬ 
tivity  to  the  presence  of  all  sizes  of  hydrometeors  but  also  resulting 
in  much  stronger  attenuation.  As  a  result  of  the  availability  of 
empirically  derived  relations,  microwave  scattering  theory,  and  the 


backscatter  coefficient 


general  absence  of  attenuation,  radar  is  at  present  a  much  more  useful 
tool  to  derive  quantitative  data. 

2.1.4  Multiple  Scattering  Correction  Factor 

For  water  and  ice  at  visible  wavelengths,  energy  absorption  is 
negligible  so  that  optical  attenuation  can  be  completely  accounted  for 
by  the  scattering  of  energy  out  of  the  laser  beam.  One-half  of  the 
total  energy  scattered  by  a  nonabsorbing,  large  sphere  is  confined  to 
a  narrow  cone  in  the  forward-scattered  direction  (i.e.,  by  diffrac¬ 
tion),  and  this  portion  of  the  scattered  energy  can  continue  to  con¬ 
tribute  to  the  illumination  of  the  target  if  the  receiver  beam  widths 
are  larger  than  the  rate  at  which  the  diffracted  component  diverges. 

The  factor  n  in  Eq.  (2.1)  accounts  for  this  diffracted  energy, 
and  ranges  between  0.5  and  1.0,  depending  on  particle  size  and  the 
lidar  beam  widths  used.  For  example,  when  looking  at  rain,  most  of 
the  diffracted  light  continues  to  propagate  in  the  forward  direction 
and  does  not  leave  the  receiver  field  of  view,  so  that  n  approaches 
0.5  (Shipley  et  al.,  1974).  However,  when  looking  at  aerosols  with 
the  lidar  (i.e.,  Rayleigh  scatterers),  there  is  no  diffraction  as 
such,  and  n  would  approach  1.0. 

2.2  Depolarization  in  the  Backscatter 
The  laser  used  in  this  study  emits  radiation  linearly  polarized 
in  the  vertical  direction.  Linear  depolarization  ratios,  or  6  values, 
are  calculated  to  characterize  the  change  in  the  state  of  polarization 


of  the  energy  backscattered  from  hydrometeors. 
The  depolarization  ratio  6  is  defined  as 
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(2.2) 


where  Pr  and  Pr  are  the  energies  returned  in  the  plane  of  polariza- 

rl  ii 

tion  orthogonal  and  parallel  to  that  of  the  transmitted  energy,  res¬ 
pectively. 

Dense  cloud  layers,  composed  only  of  cloud  droplets,  are  easily 
identified  in  lidar  returns  by  the  presence  of  near-zero  6  values  at 
the  cloud  base,  followed  by  strong  signal  increases,  evidence  for 
rapid  optical  attenuation,  and  6  values  which  gradually  increase  with 
penetration  depth.  The  latter  effect  is  caused  by  multiple  scattering 
of  the  laser  light  among  cloud  droplets. 

The  interaction,  including  scattering  and  absorption,  of  an 
electromagnetic  wave  incident  on  a  homogeneous  sphere  larger  than 
the  wavelength  can  be  described  by  using  the  Mie  solution  to  the 
Maxwell  equation.  Generally,  energy  is  backscattered  from  a  sphere 
through  two  primary  mechanisms:  surface  waves  and  paraxial  compon¬ 
ents.  For  water  spheres  and  light  in  the  visible  region,  the  results 
of  ray  tracing  methods  indicate  that  the  backscattered  energy  consists 
almost  entirely  of  rays  reflected  from  the  front  and  rear  faces  along 
the  axis  of  symmetry  (i.e.,  the  axial  rays),  while  multiple  reflected 
off-axis  rays  contribute  negligibly.  It  follows  that  the  total  axial 
intensity  varies  strongly  with  minute  variations  in  droplet  diameter 
due  to  ray  interference  effects. 

One  problem  encountered  with  the  ray  tracing  method  is  that  it 
fails  to  treat  circumferentially  backscattered  energy,  that  is,  the 


surface  wave.  This  backscattered  wave  results  when  energy  becomes 
trapped  at  the  interface  between  two  dielectrics  at  or  greater  than 
the  critical  angle  of  total  reflection.  For  water  spheres  which  are 
uniformly  irradiated,  the  axial  backscatter  component  is  small  in 
comparison  to  the  contribution  from  surface  waves  (Liou  and  Lahore, 
1974)  as  long  as  the  drop  remains  spherical,  but  a  drop  can  lose  its 
sphericity. 

If  drops  fall  at  terminal  velocity  in  the  absence  of  turbulence, 
they  are  near  perfect  spheres  if  their  equivalent  radius  (r  )  is  less 
than  140  pm.  If  their  radius  is  140  <  r  <  500  pm  they  become  oblate 
spheroids,  and  at  rm  =  500  pm  the  spheroid  develops  a  flattened  base. 
The  flattening  becomes  increasingly  pronounced  until  near  r^  =  2.0  mm 
when  a  concave  base  begins  to  appear.  As  the  drop  becomes  even 
larger,  the  concavity  deepens  until  r^  *  5  mm  the  drop  becomes  un¬ 
stable  and  breaks  up  (Pruppacher  and  Klett,  1980).  In  stratus  and 
stratocumulus  clouds,  the  average  cloud  droplet  radius  is  near  5.0  pm, 
and  the  range  of  drizzle  drops  is  between  50-500  pm.  As  a  result, 
drizzle  drops  can  become  slightly  deformed  without  turbulence  and  may 
be  a  source  of  some  depolarization. 

For  homogeneous  spherical  particles  (e.g.,  cloud  droplets  and 
small  rain  drops),  the  state  of  linear  polarization  from  sirqle  back- 
scattering  is  identical  to  that  of  the  incident  light.  However,  this 
would  not  be  true  for  nonspherical  ice  particles.  Ice  crystals,  for 
example,  cause  a  rotation  in  the  incident  plane  resulting  in  depolar¬ 
ization  due  to  multiple  internal  reflections. 

The  importance  of  irregular  ice  particle  single  scattering  cannot 
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be  overemphasized.  Some  of  the  largest  values  occur  with  large  ice 
particles  (e.g.,  graupel).  However,  this  experiment  dealt  strictly 
with  liquid  phase  clouds  of  >0°C,  so  ice  particle  scattering  will  not 
be  discussed  to  any  great  extent  here. 

Sassen  (1977b)  measured  backscattered  laser  depolarization  in  the 
laboratory  from  artificially  frozen  rain  drops  undergoing  the  solid  to 
liquid  phase  change.  When  the  drops  began  to  melt,  and  were  covered 
by  a  water  shell,  the  depolarization  values  were  relatively  high, 
almost  the  same  as  for  pure  ice.  Even  when  the  water  shell  continued 
to  increase  in  thickness  and  the  ice  core  shrank  to  the  center  of  the 
drop,  the  depolarization  values  remained  high.  The  optical  inhomo¬ 
geneity  of  the  ice  core  can  provide  opportunities  for  ray  paths 
involving  multiple  internal  reflections.  This  inhomogeneous  particle 
scattering  principle  may  have  an  analog  in  water  clouds,  as  explained 
below. 

During  the  process  of  the  condensation  of  a  clout  droplet  in  a 
maritime  environment,  the  nuclei  is  likely  to  be  a  relatively  large 
salt  particle.  Newly  formed  drops  would  have  a  salt  particle  as  the 
condensation  nuclei  at  its  core,  which  would  behave  similarly  to  the 
ice  mass  in  the  previous  example.  Even  after  the  salt  particle  is 
mostly  diluted  in  the  drop,  the  gradient  in  the  refractive  index  of 
the  drop  could  alter  the  internal  ray  paths,  perhaps  resulting  in 
increased  depolarization. 


2.3  Multiple  Scattering  in  Water  Clouds 


2.3.1  Early  Laboratory  and  Theoretical  Studies 

As  previously  discussed,  it  would  be  expected  that  no  depolariza¬ 
tion  would  result  from  backscattering  by  spherical,  homogeneous 
hydrometeors.  However,  this  has  not  been  shown  to  be  the  case.  Early 
laboratory  studies  of  artificial  water  droplet  clouds,  using  contin¬ 
uous  wave  (cw)  laser/lidar  analog  systems,  measured  linear  depolariza¬ 
tion  ratios  of  ~<5%  (Schotland  et  al.,  1971;  Liou  and  Lahore,  1974; 
and  Sassen,  1974).  It  was  known  that  the  reason  for  this  depolariza¬ 
tion  was  due  to  the  multiple  scattering  contribution  to  the  back- 
scattered  signal.  This  invalidates  the  time-independent  scattering 
assumption,  since  a  large  number  of  photon  paths  are  permissible  back 
to  the  receiver. 

Liou  and  Schotland  (1971)  used  a  computational  approach  to  cal¬ 
culate  the  6  value  from  a  volume  of  spherical  and  uniformly  distrib¬ 
uted  cloud  droplets.  They  used  two  different  cloud  models  (C4  and 
C8),  representing  fair  weather  cumulus  (rm  =  4  ym)  and  cumulus  con- 
gestus  clouds  (rm  =  8  urn),  respectively.  In  both  cases  the  drop 
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number  density  (N)  was  set  at  100  cm  and  the  vertical  probe  pene¬ 
tration  depth  was  calculated  to  a  depth  of  60  m  above  the  cloud  base. 
The  receiver  beam  width  was  also  varied.  The  maximum  6  values  were 
found  to  vary  from  3%  at  10  mrad  to  <1%  at  5  mrad  for  the  C4  model , 
and  from  4%  at  10  mrad  to  2.5%  at  5  mrad  for  the  C8  cloud  model.  In 
other  words,  6  values  increase  as  the  drop  radius  increases  for  a 
given  drop  concentration. 
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Eloranta  (1972)  used  a  small  angle  approximation  to  calculate 
the  Nth  order  contribution  to  a  lidar  return  from  a  cloud.  He 
calculated  6  as  a  function  of  penetration  depth,  and  found  maximum 
6  values  about  four  times  as  large  as  those  predicted  by  Liou  and 
Schotland  (1971),  reaching  a  maximum  of  approximately  13%  at  120  m. 

2.3.2  Early  Lidar  Measurements 

Early  laboratory  studies  were  useful  but  could  not  actually  model 
atmospheric  clouds.  Usually,  a  dense  cloud  had  to  be  created  in  a 
small  cloud  chamber.  The  cloud  had  to  be  dense  to  make  up  for  the 
short  physical  depth  of  the  chamber,  which  is  unlike  the  situation  of 
a  laser  pulse  penetrating  hundreds  of  meters  through  a  water  cloud  in 
the  atmosphere.  Pal  and  Carswell  (1973)  measured  backscattered 
depolarization  from  cumulus  using  a  ruby  laser  with  5  mrad  beam 
widths.  They  found  6  values  in  water  clouds  that  approached  15%  in 
the  first  60  m  penetration  depth,  and  depolarization  ratios  of  ap¬ 
proximately  45%  at  150  m.  These  values  were  much  higher  than  those 
predicted  by  the  early  laboratory  studies  or  theoretical  simulations, 
prompting  new  studies  to  explain  these  discrepancies. 

2.3.3  More  Recent  Studies 

Pal  and  Carswell  (1976)  again  used  cumulus  clouds  to  investigate 
multiple  scattering  in  the  atmosphere.  Since  multiple  scattering 
leads  to  depolarization  of  the  incident  pulse,  and  single  scattering 
retains  the  incident  polarization  in  the  backscatter,  they  investi¬ 
gated  the  multiple  scattered  signal  by  using  the  depolarization  as  a 
tracer.  In  other  words,  they  separated  the  contribution  of  multiple 
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scattering  from  the  return  signal  by  using  the  depolarized  component. 

Unfortunately,  for  the  sake  of  comparison  they  did  not  plot  6  as  a 

function  of  penetration  depth,  but  instead  used  t  ,  an  attenuation 

coefficient  similar  to  n»  versus  the  penetration  depth.  They  found 

that  for  a  cumulus  cloud  t  starts  at  the  cloud  base  with  a  value 

m 

close  to  zero  and  increases  with  increasing  penetration  depth. 
According  to  Pal  and  Carswell  (1976),  the  presence  at  the  cloud  base 
of  an  extremely  low  t  is  quite  characteristic  of  low  lying  cumulus 
and  would  indicate  the  presence  of  a  dilute  layer  of  single  scattering 
spherical  particles. 

Houston  and  Carswell  (1978)  investigated  atmospheric  clouds 
using  a  ruby  lidar  equipped  with  a  four-channel  receiver.  This 
receiver,  with  a  5  mrad  field  of  view,  was  capable  of  independently 
measuring  the  four  components  of  the  Stokes  vectors,  I,  Q,  U,  and  V. 
Since  6  is  a  measure  of  two  of  the  Stokes  vectors,  they  plotted  <5 
versus  the  penetration  depth  for  a  heavy  overcast  cloud  layer  with 
bases  near  700  m.  They  reported  a  maximum  6  value  of  17%  after  an 
approximately  80  m  penetration  depth.  As  far  as  the  other  components 
of  the  Stokes  vector  were  concerned,  Houston  and  Carswell  recommended 
".  .  .  that  two-channel  measurements  of  6  appear  to  be  the  most  effi¬ 
cient  way  to  utilize  lidar  for  polarization  studies." 

In  controlled  laboratory  studies,  Ryan  et  al .  (1979)  and  Sassen 
and  Liou  (1979)  investigated  <5  values  as  a  function  of  the  attenuation 
coefficient  and  optical  depth.  Ryan  et  al .  used  a  set  of  field  stops, 
placed  over  the  receiver  lens,  to  exclude  all  single  scattered  light 
from  the  receiver.  Upon  investigating  the  dependence  of  the  linear 


depolarization  ratio  versus  the  cloud  attenuation  coefficient,  they 
found  that  depolarization  increases  rapidly  with  increasing  attenua¬ 
tion,  but  levels  off  at  6  =  40%.  This  was  attributed  to  the  fact 
that  as  the  cloud  density  increases,  the  depth  of  penetration  of  the 
laser  beam  into  the  cloud  effectively  decreases.  The  very  large  6 
values  can  be  attributed  to  the  relatively  large  receiver  beam  widths 
(5  and  10  mrad),  and  to  the  blocked  primary  scattering. 

Sassen  and  Liou  (1979)  also  measured  6  as  a  function  of  cloud 

attenuation.  They  found  that  a  linear  relationship  exists  between 

increases  in  6  values  and  optical  depth.  Using  a  receiver  that  viewed 
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a  solid  angle  of  10  steradians,  they  measured  maximum  6  values  of 
~6%  when  the  extinction  coefficient  was  3.0  m_*.  Their  measurements 
were  at  e  =  1750,  as  close  to  the  backscatter  direction  as  possible  in 
the  laboratory. 

Using  a  polar  nephelometer,  they  also  investigated  the  angular 
depolarizing  patterns  for  both  water  and  ice  clouds.  Upon  investi¬ 
gating  the  relationship  between  6  and  the  scattering  angle  e  for  water 
clouds,  Sassen  and  Liou  (1979)  found  that  the  maximum  6  values  of  a 
few  percent  were  in  the  near-backward  direction  with  the  minimum  of 
-0.1%  in  the  near- forward  direction.  This  demonstrates  the  effects 
of  diffraction,  alluded  to  in  subsection  2.1.4.  For  spherical  homo- 
genoeous  particles  in  the  Mie  domain,  scattering  in  the  forward 
direction  is  governed  by  diffraction  and  grazing  reflections  which  do 
not  involve  any  depolarization. 

Carswell  and  Pal  (1980)  photographed  the  laser  backscattering 
field  from  a  water  suspension  of  monodisperse  polystyrene  spheres 


through  a  linear  polarizer.  They  found  for  sphere  sizes  in  the  Mie 
scattering  domain  that  the  backscattered  field  is  highly  polarized. 
When  the  viewing  polarizer  has  its  axes  parallel  to  the  incident 
polarization  direction,  a  cross  is  formed  with  maxima  in  the  vertical 
and  horizontal  directions.  When  the  polarizer  is  rotated  90°,  to  the 
perpendicular  direction,  the  arms  of  the  cross  rotate  45°.  Super¬ 
imposing  the  two  cases  would  create  a  complete  circle  indicating  the 
total  backscatter. 

Throughout  this  research,  the  effect  of  the  receiver  field  of 
view  was  suspected  of  causing  the  wide  range  of  values  reported.  Cai 
and  Liou  (1981)  used  numerical  methods  to  develop  a  time-dependent 
transfer  model  for  double,  triple,  and  general  multiple  backscat- 
tering.  One  of  their  more  significant  findings,  of  importance  here, 
was  the  effect  of  the  receiver  field  of  view  on  multiple  scattering 
depolarization  ratios.  They  found  that  multiple  backscattering  is 
significant  when  the  receiver  field  of  view  is  >10  mrad.  Their  cal¬ 
culations  showed  a  6  value  of  20-30%  for  a  cumulus  cloud  model  when 
the  transmitter  beam  width  and  the  receiver  field  of  view  were  10 
mrad,  and  this  was  reduced  to  <5%  when  for  1  mrad  beam  widths. 
Additionally,  6  values  decrease  to  <0.2%  for  a  larger  transmitter 
beam  width  (10  mrad)  coupled  with  a  smaller  receiver  field  of  view 
(1  mrad). 

Spinhirne  et  al .  (1983)  used  the  lidar  from  a  new  vantage  point, 
that  of  looking  down  on  cloud  tops  from  a  high  altitude  aircraft. 
Linear  depolarization  ratios  were  measured  and  found  to  increase  with 
increasing  (ice)  cloud  top  heights  and  decreasing  temperatures.  Of 


more  significance  to  this  study  was  the  amount  of  multiple  scattering 
depolarization  obtained  by  the  airborne  lidar  in  low  altitude  water 
clouds.  The  maximum  40-60%  depolarization  ratios  were  greater  than 
previously  reported  measurements  or  calculations  for  comparable 
receiver  fields  of  view  (1.5  mrad).  The  increase  was  attributed  to 
the  increased  range  factor  (15  km  versus  1  km)  for  the  airborne  obser¬ 
vations  and  the  resulting  larger  beam  footprint  diameter  at  the 
boundary  of  the  scattering  medium. 


CHAPTER  3 


STRATUS  CLOUD  PROPERTIES 


3.1  Cloud  Condensation  Nuclei 


Mixed  layers  of  stratus  and  stratocumulus  clouds  are  climato- 
logically  common  over  the  oceans  west  of  the  continents  near  the  loca¬ 
tion  of  semipermanent  high  pressure  systems.  (For  a  comparison  of 
typical  cloud  parameters,  see  Table  2.1.)  Generally,  these  layers 
can  be  divided  into  two  distinct  types  using  the  air  parcel  trajectory 
and  the  air  mass  history.  For  the  purposes  of  this  study,  the  terms 
maritime  and  continental  will  be  used  to  identify  each  distinct  type. 

The  main  difference  between  the  two  types  is  in  the  number  and 
chemical  nature  of  the  cloud  condensation  nuclei  (CCN).  Generally, 
there  is  no  latitudinal  variation  in  CCN  concentration,  with  contin¬ 
ental  air  masses  having  more  CCN  than  maritime  air  masses.  The  CCN 
concentration  normally  decreases  with  height  over  the  continents.  In 
contrast,  maritime  CCN  concentrations  may  remain  fairly  constant  with 
height  and  sometimes  may  even  increase  above  the  marine  inversion 
(Pruppacher  and  Klett,  1980). 

The  different  chemical  nature  of  the  CCN  leads  to  a  difference 
in  cloud  droplet  spectra  n(r)  (i.e.,  number  per  radius  interval)  and 
liquid  water  content,  LWC.  Noonkester  (1984)  measured  n(r)  in  marine 
stratus  layers  110-150  km  southwest  of  San  Diego  in  May  and  August 
1981,  and  found  that  the  droplet  spectra  on  the  days  influenced  by 
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continental  air  contained  many  more  small  droplets  than  on  maritime 
days,  a  fact  also  found  by  Goodman  (1977). 

Ryan  et  al .  (1972)  used  an  aircraft-mounted  particle  spectrometer 
probe  to  measure  cloud  droplet  distributions  in  California  stratus 
30-60  km  west  of  Santa  Cruz.  As  shown  in  Figure  3.1,  although  the 

.3 

total  droplet  concentration  (N,  cm  )  remains  fairly  constant,  the 
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liquid  water  content  (LWC,  gm  )  increases  and  the  drop  distribution 
becomes  broader  while  ascending  from  the  base  to  the  top  of  the  cloud. 
This  cloud  was  described  as  a  ".  .  .  shallow  non-precipitation  cloud 
whose  nucleus  spectrum  is  typically  maritime." 

Goodman  (1977)  collected  data  from  an  instrumented,  250  m  high 
tower  located  atop  a  250  m  hill  in  the  center  of  San  Francisco.  The 
coastal  stratus  was  advected  onshore  by  westerly  flow,  which  is 
similar  to  the  situation  at  the  test  site  used  in  this  study.  In 
summary,  it  was  found  that  the  drop  size  distribution  was  a  function 
of  the  vertical  distance  from  the  inversion  base,  with  the  mean 
diameter  increasing  with  height  (see  Figure  3.2).  Furthermore,  the 
drop  size  distribution  was  found  to  correlate  highly  with  the  droplet 
concentration.  Those  cases  with  a  high  concentration  had  a  narrow 
spectrum  with  a  well-defined  peak  and  a  relatively  small  mean  diam¬ 
eter.  That  would  be  expected,  since  the  available  water  vapor  must 
be  distributed  over  more  nuclei.  The  development  of  drizzle-sized 
drops  was  observed  in  all  cases  in  the  uppermost  part  of  the  stratus, 
and  later  in  the  life  cycle  of  the  cloud,  precipitation-sized  drops 
occasionally  reached  the  ground,  as  they  did  in  this  study. 

Goodman  concluded  that  maritime  (westerly)  trajectories  resulted 
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Fig.  3.1  Cloud  particle  distributions  for  a  gradual  ascent  through 
shallow  nonprecipitating  maritime  stratus.  Z  represents 
height  above  cloud  base  (after  Ryan  et  al.,  1972). 
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in  low  droplet  concentrations  (89  cm  1,  large  mean  diameters,  and 
broad  drop  size  distributions.  However,  continental  (northwest 
through  north)  trajectories  have  higher  droplet  concentrations  (265 
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cm  ),  smaller  mean  diameters,  and  narrow  drop  size  distributions  with 
sharper  peaks. 


3.2  Cloud  Water  Content 

According  to  Pruppacher  and  Klett  (1980),  the  distribution  of 
the  cloud  water  content  has  four  characteristic  features.  Two  of 
these  characteristics  which  are  pertinent  to  this  discussion  are: 

1)  the  cloud  water  content  increases  with  height  above  cloud  base, 
assumes  a  maximum  somewhere  in  the  upper  half  of  the  cloud,  and  then 
decreases  again  toward  the  cloud  top;  and  2)  the  distribution  of  the 
water  content  parallels  the  distribution  of  the  drop  size  rather  than 
the  drop  concentration. 

The  mean  LWC  measured  by  Ryan  et  al .  (1972)  was  0.20-0.35  gm“3. 
The  liquid  water  content  found  by  Goodman  (1977)  was  calculated  dir¬ 
ectly  from  the  drop  size  distributions,  and  is  shown  in  Figure  3.2. 


CHAPTER  4 


THE  FIELD  EXPERIMENT 

4.1  Location 

The  field  site  was  located  at  San  Juan  Capistrano,  California, 
at  the  TRW  Corporation  test  site  185  m  above  mean  sea  level  (MSL). 

The  experiment  dates  were  June  through  August  1982.  The  atmospheric 
soundings  used  in  this  research  were  taken  at  Montgomery  Field  in 
San  Diego,  30  km  southeast.  These  locations  are  shown  in  Figure  4.1. 

4.2  Equipment 


4.2.1  Lidar  System 

The  lidar  system  consists  of  a  transmitter,  a  pulsed  ruby  laser 
with  collimator,  and  the  receiver,  a  25  cm  diameter  telescope.  The 
transmitter  and  receiver  are  monostatic,  located  side  by  side  with 
parallel  axes  so  that  the  scattered  signal  is  -180°  with  respect  to 
the  direction  of  the  incident  beam.  The  entire  system  is  mounted  on 
a  table  supported  by  a  searchlight  mount  with  manual  azimuth  and 
elevation  angle  control.  The  mount  is  on  a  track  inside  a  van,  so 
that  the  lidar  system  can  be  positioned  inside  or  outside  of  the  van. 

The  receiver  design  incorporates  a  polarizing  prism  and  dual 
photodetectors  to  simultaneously  measure  the  components  of  the 
returned  laser  energy.  Since  multiple  scattering  strongly  depends  on 
the  receiver  field  of  view,  the  receiver  beam  width  can  be  rapidly 


San  Diego 


Fig.  4.1  Map  of  southern  California  showing  the  lidar  field  site 
and  local  rawinsonde  release  locations. 


The  lidar  data  collection  and  analysis  were  accomplished  using  a 


dual-channel  digital  oscilloscope  with  a  magnetic  disk  recorder  in 
combination  with  a  laboratory-based  minicomputer.  Lidar  returns  were 
obtained  at  0.25  to  5  min  intervals,  depending  on  what  particular 
experiment  was  being  performed,  by  simultaneously  digitizing  the 
parallel  and  perpendicular  polarized  signals  at  a  rate  of  50  ns  per 
point,  yeilding  a  spatial  resolution  of  7.5  m  in  range.  The  linear 
depolarization  ratios  have  since  been  calculated  for  each  pair  of  data 
points  after  compensating  for  differences  in  the  receiver  gain 
channels  determined  through  a  calibration  at  the  field  site. 

4.3  Synoptic  Situation 

Discussed  in  this  and  the  following  sections  will  be  some  general 
climatology  and  the  synoptic  situation  for  southern  California  during 
the  experiment.  Included  along  with  these  discussions  are  the 
National  Weather  Service  (NWS)  1200  GMT  surface  weather  maps,  geo¬ 
stationary  satellite  photographs,  and  sounding  data  from  the 
Montgomery  Field  site. 

Southern  California  generally  experiences  a  Mediterranean  type 
climate  characterized  by  warm  to  hot  summers  and  mild  winters.  Cold 
fronts  rarely  penetrate  this  area,  with  most  of  the  precipitation 
resulting  from  moist  Pacific  air  being  forced  up  the  slopes  of  the 
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Table  4.1  Ruby  lidar  specifications. 


Transmitter 


Wavelength 

Energy/pulse 

Pulse  repetition  frequency 
Pulse  duration 
Spatial  resolution 
Polarization 
Beam  width 

Receiver 


Aperture  diameter 

Polarization 

Detector 

Beam  width 

Positioning 

Data  Handling 

Digitizer  resolution 
Sample  width 
Samples/return 
Sampling  rate 
Real-time  display 
Recording 


.694  pm 

1.5  J  maximum 
0.2  Hz  maximum 
3x10- 8  s 

4.5  m 
Vertical 
1x10" 3  rad 


0.25  m 

Vertical  and  horizontal 
Dual  photomulitplier  tubes 
1-3x10" 3  rad,  variable 
Steerable  mount 


8  bit  (2  channel) 

7.5  m 

2,000  per  channel 
20  pulses  per  minute 
Cathode  ray  tube 
Magnetic  disk 


Coastal  Range.  During  the  summer,  as  the  subtropical  high  pressure 
area  moves  north,  extremely  stable  air  moves  into  the  region  from  the 
west  and  northwest,  causing  widespread  stratus,  stratocumulus,  and 
fog  situations. 

4.3.1  Daily  Weather  Discussion 

The  NWS  1200  GMT  surface  chart  for  10  June  1982  (see  Figure  4.2) 
shows  a  low  pressure  center  located  near  Las  Vegas,  Nevada  with  an 
inverted  trough  extending  along  the  west  side  of  the  Coastal  Range 
in  a  SE-NW  direction.  The  500  mb  chart  indicates  convergence  over 
the  San  Diego  area  with  a  moderate  southwesterly  flow.  During  the 
period  between  10-14  June,  the  surface  chart  showed  a  continued 
northerly  migration  of  the  subtropical  high  with  most  of  California 
dominated  by  a  weak  NW-N  (continental)  flow.  On  14  June  a  weak  cold 
front  was  located  between  Los  Angeles  and  San  Diego  with  the  500  mb 
flow  from  the  northwest.  On  15  June  there  was  a  return  to  the 
pattern  established  on  10  June,  with  an  inverted  trough  along  the 
mountains  and  a  low  pressure  to  the  east. 


4.3.2  Satellite  Imagery 

The  satellite  photograph  of  Figure  4.3  for  10  June  1982  shows  a 
widespread  area  of  stratocumulus  extending  from  the  Coastal  Range  to 
hundreds  of  kilometers  into  the  Pacific.  In  the  immediate  vicinity 
of  San  Diego,  the  cloud  tops  appear  to  be  more  stratiform,  becoming 
more  cumuliform  (i.e.,  cellular)  to  the  west.  Other  satellite  photo¬ 
graphs  from  14-15  June  are  remarkably  similar  to  Figure  4.3,  al¬ 
though  the  smoother  appearing  tops  extend  further  to  the  west. 


4.3.3  Atmospheric  Soundings 

The  sounding  from  10  June,  selected  as  typical  for  the  observa¬ 
tion  period,  indicates  two  different  inversions  above  the  ground  level 
(AGL)  of  the  lidar  field  site.  The  first  is  a  weak  frontal  inversion 
between  0.8- 1.0  km,  with  a  subsidence  inversion  present  between  1.0- 
2.0  km.  The  lifting  condensation  level  (indicative  of  the  cloud  base) 
is  at  0.4  km.  Above  1.0  km  the  sounding  becomes  extremely  dry,  indi¬ 
cating  the  upper  limit  of  the  cloud  layer.  Winds  throughout  the  lower 
1.0  km  are  light  and  westerly  and  become  moderate  southwesterl ies 


above  the  inversion  at  2.0  km  (see  Figure  4.4). 


0 

TEMPERATURE  (C) 


Fig.  4. 


4  Atmospheric  sounding  data  for  the  morning  of  10  June  1982, 
from  Montgomery  Field.  The  height  above  ground  level  (AGL) 
for  the  remote  sensor  site  (185  m  MSL)  is  shown. 


CHAPTER  5 


EXPERIMENTAL  RESULTS 

5.1  Temporal  Variations  During  Stratus  Breakup 
The  Height-Time  (HTI)  display  of  Figure  5.1  shows  a  layer  of 
stratus  undergoing  lifing  and  dissipation  on  10  June  1982.  The 
ordinate  is  divided  into  hundreds  of  meters  above  ground  level  (185 
MSL),  and  the  abscissa  is  in  local  time  (10  minute  increments).  The 
apparent  cloud  bases  and  tops  are  indicated  by  a  solid  line  and  a 
dashed  line,  respectively.  Areas  of  strong  aerosol  returns,  appar¬ 
ently  left  over  from  cloud  dissipation,  are  indicated  by  short  dashed 
lines.  Linear  depolarization  ratio  intervals,  in  percent,  are  shaded 
and  identified  in  the  key,  and  range  from  near  0  to  22.4%.  All  of  the 
data  for  this  HTI  were  obtained  at  a  90°  elevation  angle,  using  a 
1  mrad  receiver  field  of  view. 

Before  continuing  with  the  discussion  of  the  HTI,  a  discussion  of 
how  the  cloud  base  was  determined  is  warranted.  When  referring  to 
the  cloud  base  from  the  atmospheric  sounding  data,  the  height  of  the 
lifting  condensation  level  (LCL)  was  used.  The  LCL  is  the  height  at 
which  a  parcel  of  air  becomes  saturated  when  lifted  dry  adiabatical ly. 
When  referring  to  the  lidar  cloud  base,  a  determination  has  to  be 
made  as  to  the  exact  definition  of  a  cloud  base.  Since  the  lidar 
signals  are  affected  by  the  air  molecules  and  the  cloud  condensation 
nucleus  activation  zone  beneath  a  cloud  layer,  the  backscattering 


versus  time  display  of  lidar  6  values  during  the  dissipation  of  a  stratus 
t  the  indicated  times  on  10  June  1982. 
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signals  typically  begin  to  increase  below  the  level  corresponding  to 
100%  relative  humidity,  particularly  in  the  maritime  environment. 

For  this  experiment,  the  cloud  base  was  defined  as  the  height  at  which 
the  first  of  three  consecutive  10%  parallel-polarized  energy  increases 
were  measured  in  the  digitized  signals. 

5.1.1  Visual  Cloud  Observations 

Since  optical  surface  observations  were  unavailable  for  this 
specific  location,  it  will  be  helpful  to  describe  the  visual  observa¬ 
tions  of  the  stratus  layer  during  dissipation. 

At  1130  LT  (Local  Time),  the  lidar  operator  reported  a  "ridged” 
stratus  overcast  with  intermittent  drizzle  occurring.  Thirty  minutes 
later  a  few  blue  patches  appeared  before  returning  to  the  overcast 
condition.  At  1240  LT  the  operator  reported  that  the  sun  was  visible 
with  7/10  cloud  coverage,  and  at  1310  LT,  4/10  "cumulus-type"  clouds 
decreasing  to  2/10  coverage  at  1330  LT  were  reported. 

5.1.2  Lidar  Findings 

At  the  beginning  of  the  time  period,  between  1010-1030  LT,  the 
stratus  cloud  base  was  indicated  at  approximately  375  m,  with  ragged 
tops  extending  to  approximately  600  m  above  ground  level  (AGL).  The 
cloud  top  heights  shown  in  Figure  5.1  may  not  represent  the  actual 
top  of  the  layer,  however,  due  to  the  attenuation  of  the  laser  pulse. 
However,  the  maximum  cloud  heights  roughly  agree  with  the  sounding 
data  for  this  occasion. 

A  good  agreement  exists  between  the  operator's  observations, 
the  lidar  measurements  of  the  stratus  layer,  and  the  cloud  heights 
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derived  from  the  atmospheric  sounding  of  Figure  4.4,  which  indicated 
a  possible  cloud  layer  between  0. 4-1.0  km  AGL.  Throughout  the  period, 
the  stratus  layer  lifted  with  intermittent  drizzle  occurring,  espe¬ 
cially  at  1130  LT.  At  that  time,  the  lidar  indicated  a  depolarization 
value  of  <2.4%  below  the  cloud  base.  A  value  this  low  could  only  be 
due  to  spherical,  homogeneous  droplets  at  relatively  low  concentra¬ 
tions.  In  other  words,  light  drizzle  was  indicated  below  the  cloud 
base.  The  effect  of  multiple  scattering  or  depolarization  under  these 
conditions  is  minimal  due  to  the  1  mrad  receiver  field  of  view.  The 
period  after  1240  LT  indicated  a  break-up  of  the  stratus  layer  with 
holes  between  the  cloud  elements.  The  periods  near  1310  LT  and  1335 
LT  indicated  an  enhanced  aerosol  return,  probably  attributable  to  the 
dissipation  of  the  stratus  clouds.  This  area  had  a  depolarization 
value  less  than  the  surrounding  cloud  elements,  but  greater  than  the 
weak  aerosol  returns  below  the  cloud  base  (generally  between  2. 5-7. 5%) 
which  have  been  omitted  from  the  figure. 

Above  the  cloud  base,  the  value  of  the  depolarization  ratio 
typically  increases  rapidly  with  increasing  penetration  depth.  Higher 
in  the  cloud,  however,  the  6  value  variations  are  clearly  more  erratic 
with  position  and  time.  This  is  due  to  multiple  scattering  effects, 
coupled  with  cloud  microphysical  content  variations.  Some  of  the 
microphysical  variations  to  be  expected  in  this  cloud  layer  during 
dissipation  can  be  seen  by  referring  to  Figure  3.2  (after  Goodman, 
1977). 

Figure  3.2  shows  the  properties  of  a  stratus  layer  covering  the 
period  from  formation  to  dissipation.  It  illustrates  that  the  liquid 


water  content  and  mean  drop  size  distribution  have  increasing  values 
with  height  and  time  until  reaching  the  dissipation  stage.  Then  the 
values  for  both  LWC  and  drop  size  decrease  below  the  maxima  attained 
in  the  middle  stages.  The  droplet  concentration  has  no  noticeable 
trend  except  for  a  pronounced  decrease  in  the  upper  portion  of  the 
cloud  layer. 

The  relation  of  these  general  cloud  properties  to  the  depolar¬ 
ization  data  is  unclear,  but  early  in  the  period  shown  in  Figure  5.1, 
when  the  apparent  cloud  depth  was  large,  the  range  of  6  values  was 
small,  only  2.5%  to  10%.  However,  as  the  cloud  layer  appeared  to 
thin,  the  range  of  6  values  increased  from  2.5%  to  17.5%,  with  the 
largest  area  of  high  6  values  occurring  at  1300  LT  when  6  reached 
values  between  17.5%  and  22.4%.  Cloud  composition  changes  and  the 
increasing  distance  to  the  target  may  both  have  contributed  to  the 
increasing  6  values  from  multiple  scattering  with  time  and  height. 

5.2  Receiver  Beam  Width  Effects 

Throughout  this  research,  the  effect  of  the  lidar  receiver  beam 
width  has  been  identified  as  a  major  contributing  factor  to  depolar¬ 
ization  ratios  due  to  multiple  scattering  in  water  clouds.  As 
reported  earlier,  water  clouds  were  probed  with  lasers,  modeled, 
numerically  simulated,  and  found  to  have  6  ratios  as  low  as  a  few 
percent  to  as  high  as  almost  50%.  The  receiver  field  of  view  was 

ususally  different  from  one  experiment  to  the  next,  ranging  from 
-2  -4 

10  rad  to  10  rad.  Even  though  different  results  for  the  same 
beam  widths  have  been  reported,  no  study  has  systematically  varied 
the  receiver  beam  width  and  used  the  lidar  in  the  atmosphere  to 
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measure  the  backscattered  depolarization  from  water  clouds.  The 
following  nine  figures  represent  the  first  time  that  this  has  been 
accomplished  (Figures  5.2  -  5.10). 

These  figures,  from  14  June,  will  be  divided  into  groups  accord¬ 
ing  to  the  lidar  elevation  angle.  They  were  selected  as  typical 
examples  from  a  sample  of  approximately  one  hundred  lidar  returns  for 
this  day.  There  are  four  figures  at  90°,  three  at  60°,  and  one  each 
at  50°  and  30°.  Most  of  these  figures  depict  the  linear  depolariza¬ 
tion  ratio  (in  %>)  as  a  function  of  vertical  penetration  depth  (m) 
above  the  cloud  base,  and  compare  1  mrad  and  3  mrad  receiver  fields 
of  view.  Beneath  each  6  value  plot  is  the  corresponding  graph  of  the 
relative  returned  energy  in  the  parallel  polarization  plane,  given  in 
arbitrary  units.  Unlike  the  HTI  data  display  of  Figure  5.1,  which 
has  the  appearance  of  an  actual  cloud  layer,  these  figures  are  der¬ 
ived  from  single  lidar  pulses  and  resemble  the  A-scope  displays 
available  in  real  time  from  most  modern  radar  or  lidar  systems. 

In  Figure  5.2,  instead  of  using  the  penetration  depth  above  cloud 
base,  the  depolarization  ratios  and  relative  returned  energies  are 
shown  as  a  function  of  height  above  ground  level.  These  have  been 
added  to  illustrate  the  apparence  of  a  lidar  return  through  the  atmo¬ 
sphere  from  the  surface  to  approximately  600  m.  The  specific  numer¬ 
ical  values  from  this  figure  are  not  important  at  this  time,  but  an 
understanding  of  what  the  lidar  senses  and  how  that  will  be  displayed 
is  important.  This  figure  shows  a  low  altitude  aerosol  plus  molec¬ 
ular  return  and  a  stratus  cloud  return  beginning  at  500  m  AGL  for 
both  3  and  1  mrad  receiver  fields  of  view.  The  rapid  increase  in 


HEIGHT  AGL  CM)  HEIGHT  AGL  CM) 


DEPOLARIZATION  RATIO 


RELATIVE  ENERGT 


Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  height  above  ground  level  for  two  consecutive  lidar 
pulses  using  1  and  3  mrad  (MR)  receiver  beam  widths.  The 
lidar  elevation  angle  was  90°.  The  data  were  collected 
at  1036  LT  on  14  June  1982. 
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Fig.  5.4  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  two  con 
secutive  lidar  pulses  using  1  and  3  mrad  (MR)  receiver 
beam  widths.  The  data  were  collected  at  a  90°  elevation 
angle  at  1259  LT  on  14  June  1982  from  a  cloud  with  a  480  m 
cloud  base. 
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Fig.  5.6  Linear  depolarization  ratio  (in  %)  and  relative  returned 

energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  two  con 
secutive  lidar  pulses  using  1  and  3  mrad  (MR)  receiver 
beam  widths.  The  data  were  collected  at  a  60°  elevation 
angle  at  1124  LT  on  14  June  1982  from  a  cloud  with  a  450  m 
cloud  base. 
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Fig.  5.7  Linear  depolarization  ratio  (inland  relative  returned 

energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  two  con 
secutive  lidar  pulses  using  1  and  3  mrad  (MR)  receiver 
beam  widths.  The  data  were  collected  at  a  60°  elevation 
angle  at  1217  LT  on  14  June  1982  from  a  cloud  with  a  410  m 
cloud  base. 
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Fig.  5.8  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  two  con 
secutive  lidar  pulses  using  1  and  3  mrad  (MR)  receiver 
beam  widths.  The  data  were  collected  at  a  60°  elevation 
angle  at  1448  LT  on  14  June  1982  from  a  cloud  with  a  630  m 
cloud  base. 
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Fig.  5.10  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  two  con 
secutive  lidar  pulses  using  1  and  3  mrad  (MR)  receiver 
beam  widths.  The  data  were  collected  at  a  30°  elevation 
angle  at  1418  LT  on  14  June  1982.  The  penetration  depth 
scale  begins  at  540  m  AGL. 


returned  energy  below  120  m,  the  point  at  which  the  transmitter/ 
receiver  beam  widths  overlap  with  this  lidar  (see  Section  2.1.2),  is 
indicative  of  multiple  scattering  among  air  molecules  and  aerosols. 
Above  this  level  the  returned  energy  tends  to  decrease  due  to  the 
inverse  range-squared  law,  until  a  gradual  and  then  sudden  increase 
marks  the  cloud  nuclei  activation  zone  and  the  strong  cloud  base 
return.  Another  notable  feature  pertains  to  the  depolarization  ratio. 
It  is  significantly  higher  for  the  3  mrad  as  opposed  to  the  1  mrad 
field  of  view,  as  will  be  considered  in  the  next  section. 

5.2.1  Observations  at  a  90°  Elevation  Angle 

To  begin  the  discussion  of  the  variable  beam  width  data  collected 
on  14  June  1982,  Figure  5.3,  from  1234  LT,  shows  a  maximum  depolariza¬ 
tion  ratio  about  three  times  as  large  for  3  mrad  (31%  versus  11%)  as 
for  1  mrad.  Figure  5.4,  obtained  25  minutes  later,  also  shows  a 
trebled  <5  value  for  the  3  versus  1  mrad  receiver  beam  widths  (5.5% 
versus  16%).  This  was  the  rule  for  most  of  the  approximately  one 
hundred  cases  reviewed  for  this  day.  When  comparing  the  3  versus  2 
mrad  data  in  Figure  5.5  (1259  LT),  the  maximum  6  values  are  nearly 
doubled.  The  maximum  relative  returned  energy  varies  roughly  between 
500-600  (arbitrary  units),  with  the  3  mrad  returns  generally  having 
a  higher  peak  value  than  1  mrad.  In  Figure  5.5,  a  double  cloud 
layer  may  have  been  encountered,  since  the  relative  energy  returned 
increased  at  20  m  and  then  decreased  before  increasing  again  at  30  m. 

Maximum  depolarization  values  at  a  90°  elevation  angle  on  this 


day  ranged  from  16%  to  slightly  over  30%  at  3  mrad,  and  between  6% 
and  15%  at  1  mrad.  It  is  apparent  that  these  data  would  fall  somewhere 
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in  the  middle  of  the  '5%  to  50%  range  reported  earlier  in  laboratory. 


theoretical,  and  atmospheric  measurements. 

5.2.2  Observations  at  a  60°  Elevation  Angle 

Figures  5.6,  5.7,  and  5.8  are  from  1124  LT,  1217  LT,  and  1448  LT, 
respectively.  For  each  of  these  figures  the  slant  range  of  the  pulse 
into  the  cloud  was  converted  to  a  vertical  penetration  depth  above 
the  cloud  base  scale  through  the  use  of  the  sine  law.  As  at  the  90° 
elevation  angle,  maximum  6  values  are  two  to  three  times  greater  for 
the  3  mrad  as  opposed  to  the  1  mrad  observations.  One  characteristic 
at  the  60°  elevation  angle  is  that  6  reaches  its  maximum  value  lower 
in  the  cloud  for  1  mrad  than  for  3  mrad,  and  usually  maintains  this 
value  until  the  signal  drops  off.  This  was  not  the  case  at  90°,  where 
the  6  profiles  for  both  1  mrad  and  3  mrad  were  more  bell  shaped.  At 
a  lower  elevation  angle,  the  laser  beam  takes  a  longer  path  to  reach 
the  same  height  in  the  cloud  as  a  result  of  its  slant  path.  Since 
the  composition  of  the  cloud  would  be  expected  to  change  above  the 
cloud  base,  lidar  returns  varying  as  a  function  of  elevation  angle 
would  also  be  expected. 

5.2.3  Observations  at  50°  and  30°  Elevation 
Angles 

The  example  of  a  50°  elevation  angle  return  in  Figure  5.9  again 
has  an  almost  threefold  increase  for  the  depolarization  ratio  between 
3  and  1  mrad.  The  shape  of  the  1  mrad  6  profile  has  a  definite  peak, 
as  in  the  90°  returns,  but  the  6  values  do  not  continue  to  decrease 
as  the  signal  decreases. 
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The  unusual  set  of  lidar  returns  given  in  Figure  5.10  is  for  30° 
elevation  angle  observations  collected  at  1418  LT.  Figure  5.10  has 
approximately  the  same  <5  values  for  both  3  and  1  mrad  up  to  -100  m 
penetration  depth,  which  are  due  to  the  presence  of  a  weak  subcloud 
layer  probably  containing  cumulofractus  clouds.  The  depolarization 
ratio  and  relative  energy  profiles  both  indicate  a  much  stronger 
scattering  cloud  layer.  Unfortunately,  there  were  no  other  observa¬ 
tions  performed  at  30°  at  this  time. 

5.3  Elevation  Angle  Effects 

In  the  previous  sections,  elevation  angle  effects  naturally 
entered  into  the  description  of  the  receiver  beam  width  effects. 
However,  the  time  between  the  observations  at  various  elevation 
angles  was  often  quite  large,  making  it  difficult  for  intercompari¬ 
sons.  For  the  next  three  figures,  the  time  between  shots  varies  only 
by  minutes,  so  the  elevation  angle  effects  can  be  more  closely 
compared. 

Figure  5.11  compares  90°,  60°,  and  30°  elevation  angle  observa¬ 
tions  collected  with  a  3  mrad  receiver  beam  width  at  1427  LT,  1429 
LT,  and  1431  LT,  respectively.  The  patterns  discussed  in  the  pre¬ 
vious  sections  still  hold  true  for  these  figures.  For  example,  maxi¬ 
mum  depolarization  values  and  relative  returned  energy  peaks  tend  to 
be  reached  lower  in  the  cloud  as  the  angle  is  decreased. 

Figure  5.12  is  given  to  assess  the  temporal  variability  of  the 
data  shown  in  Figure  5.11.  These  data  were  collected  one  minute 
later  at  each  angle.  Both  figures  are  very  similar,  showing  the  con¬ 
servative  properties  of  the  stratus  layer  and  the  lidar  system  over  a 
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11  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  3  mrad 
observations  taken  2  min  apart  at  the  indicated  elevation 
angles.  Average  cloud  base  was  at  545  m  AGL. 
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Fig.  5.12  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  vertical  penetration  depth  above  cloud  base  for  3  mrad 
observations  taken  2  min  apart  at  the  indicated  elevation 
angles.  Average  cloud  base  was  at  545  m  AGL.  The  data 
were  collected  one  minute  later  at  each  elevation  angle 
than  in  Figure  5.11. 


short  time  period. 


Figure  5.13  presents  the  data  in  Figure  5.12  as  a  function  of 
slant  path  range  into  the  cloud.  It  can  be  seen  that  both  the  maximum 
6  values  and  the  relative  energy  peaks  more  closely  coincide  when 
using  the  slant  path  distance  as  opposed  to  the  vertical  penetration 
depth  above  cloud  base.  However,  the  dissimilar  6  values  in  the 
cloud  between  the  90°  or  60°  and  30°  data  point  to  the  presence  of 
cloud  microstructure  changes  as  a  function  of  distance  above  the  cloud 
base. 


5.4  Beam  Misalignment  Effects 

This  section  of  the  experimental  results  deals  with  the  deliber¬ 
ate  misalignment  of  the  lidar  transmitter  and  receiver  axes.  As 
mentioned  in  Section  2.1.2,  the  transmitter  and  receiver  axes  must  be 
aligned  in  parallel  by  sighting  a  fixed  target  through  the  telescope 
while  steering  the  laser  pulse  with  the  collimator.  To  test  the 
effect  of  misalignments,  the  laser  beam  alignment  was  deliberately 
moved  off-axis  with  respect  to  the  transmitter.  Figure  5.14  shows 
the  effect  of  this  misalignment,  and  was  taken  at  matched  1  mrad  beam 
widths  with  a  45°  elevation  angle.  The  beam  axes  were  made  to  diverge 
in  steps  starting  from  an  angle  of  0.0  (the  aligned  case)  to  7.5  mrad 
off  alignment  in  the  vertical  direction. 

For  the  0.0  mrad  return  of  Figure  5.14,  a  cloud  layer  is  indi¬ 
cated  between  320  and  400  m  with  a  maximum  6  value  of  about  5%.  As 
the  off-axis  misalignment  increases  to  0.9  mrad,  the  returned  laser 
energy  begins  to  decrease  and  a  small  increase  in  the  6  value  is 
noted.  However,  the  next  step  (2.0  mrad  off-axis)  and  each  successive 
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Fig.  5.13  Linear  depolarization  ratio  (in  %)  and  relative  returned 
energy  (in  arbitrary  units)  profiles  plotted  as  a  function 
of  slant  path  range  into  the  cloud  for  3  mrad  observations 
taken  2  min  apart  at  the  indicated  elevation  angles. 
Average  cloud  base  was  at  545  m  A6L.  The  data  were  col¬ 
lected  one  minute  later  at  each  elevation  angle  than  in 
Figure  5.11. 
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Fig.  5.14  Effects  of  transmitter/receiver  beam  misalignment  on 

returns  from  a  stratus  cloud  layer  with  a  cloud  base  at 
320  m  above  ground  level.  Numbers  in  the  keys  refer  to 
receiver  pointing  angle  misalignment,  in  mrad.  The  ele 
vation  angle  was  45°. 


step  thereafter  shows  a  dramatic  change.  As  the  angle  of  misalign¬ 
ment  increases,  the  transmitter  and  receiver  beams  increasingly  fail 
to  overlap,  and  the  detection  of  multiple  scattered  energy  is  in¬ 
creasingly  favored  over  single  scattering.  Eventually,  the  return 
is  due  only  to  multiple  scattering. 

Regarding  the  relative  energy  portion  of  the  figure,  the  energy 
peak  returned  decreases  for  each  successive  off-axis  plot  and  moves 
further  up  in  the  atmosphere.  It  appears  that  as  the  beams  move 
further  apart  progressively  higher  orders  of  multiple  scattering  are 
detected  by  the  receiver.  These  observations  illustrate  the  time- 
dependent  nature  of  the  multiple  scattering  process,  since  the  cloud 
layer  appears  to  move  to  higher  altitudes  with  increasing  misalign- 


CHAPTER  6 


CONCLUSIONS 

The  overall  purpose  of  this  study  was  to  investigate  the  effects 
of  multiple  scattering  in  water  clouds,  and  more  specifically,  to 
test  some  hypotheses  about  the  effects  of  receiver  beam  widths  on 
depolarization  ratios  in  stratus  clouds.  This  study  represents  a  step 
in  a  progression  of  studies  aimed  toward  the  goal  of  better  under¬ 
standing  the  structure  of  water  clouds.  The  lidar  used  in  this  study, 
a  relatively  new  device,  is  an  instrument  which  will  undoubtedly  find 
increasing  applications  in  atmospheric  research.  The  entire  instru¬ 
ment  and  related  data  processing  equipment  fit  into  the  back  of  a  van 
for  easy  transport  and  can  be  on  line  and  taking  data  in  a  single 
hour. 

Cloud  droplet  induced  depolarization  from  multiple  scattering, 
first  thought  to  jeopardize  cloud  phase  discrimination  with  lidar, 
actually  contains  information  on  the  microphysical  characteristics  of 
the  cloud.  It  should  be  emphasized  that  by  choosing  a  sufficiently 
small  receiver  beam  width,  the  effects  on  the  depolarization  ratio  of 
multiple  scattering  can  be  minimized.  When  choosing  a  small  receiver 
beam  width,  however,  great  care  must  be  taken  in  aligning  the  system. 
For  receiver  beam  widths  on  the  order  of  1  mrad,  a  relatively  small 
misalignment  in  the  transmitter/ receiver  axes  can  eliminate  all  single 
scattering  and  leave  only  the  multiple  scattered  component. 
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In  the  near  future  the  FIRE  (First  International  Satellite  Cloud 
Climatology  Project  Regional  Experiment)  project  will  investigate 
marine  stratus  and  stratocumulus,  and  cirrus  clouds  to  help  assess 
their  effects  on  the  earth's  energy  budget  and  climate  system.  Laser 
systems,  used  in  conjunction  with  a  wide  range  of  other  research 
equipment,  are  scheduled  to  be  a  major  component  of  the  FIRE  experi¬ 
ments. 

During  the  study  described  here,  lidar  measurements  have  been 
reported  which  suggest  that  stratus  cloud  microphysical  variations 
may  produce  characteristic  changes  in  the  depolarization  ratio.  It 
is  recommended  that,  in  order  to  assess  the  full  information  content 
of  polarization  lidar  observations  in  stratus  clouds,  a  combined 
program  of  lidar  measurements,  in  situ  observations  of  the  internal 
cloud  structure,  and  new  multiple  scattering  theoretical  simulations 
for  inhomogeneous  clouds  be  pursued. 
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